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Abstract 

We present an analysis of the preliminary experimental data of direct J/b 
production in e + e _ process at low energies. We find that the color-octet con¬ 
tributions are crucially important to the cross section at this energy region, 
and their inclusion produces a good description of the data. By fitting to 
the data, we extract the individual values of two color-octet matrix elements: 
(ot{ x S Q )) » 1.1 x 10 ~ 2 GeV 3 , (ot( 3 P 0 )) /m 2 c « 7.4 x 10~ 3 GeV 3 . We discuss 
the allowed range of the two matrix elements constrained by the theoreti¬ 
cal uncertainties. We find that ( 1 *S'o)) is poorly determined because it is 

sensitive to the variation of the choice of m c , a s and (of( 3 Si)^. However 

((b^( 3 Po)) /m( is quite stable (about (6 — 9) x 10 -3 GeH 3 ) when the param¬ 
eters vary in reasonable ranges. The uncertainties due to large experimental 
errors are also discussed. 

PACS number(s): 13.60.Le, 13.20.Gd 
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In recent years, NRQCD factorization formalism [|T| has made impressive progress in 
the study of heavy quarkonium production j2j. In this approach, the production process is 
factorized into short and long distance parts, while the latter is associated with the non- 
perturbative matrix elements of four fermion operators. The cross section for the inclusive 
production of a quarkonium state H can be expressed as a sum of products having the form 

da(A + B^H + X) = J2 da(A + B^ cc[n} + X ) (O*) . (1) 

n 

In the above, da represent the short distance coefficients associated with the perturbative 
subprocesses in which a cc pair is produced in a configuration denoted by n (angular mo¬ 
mentum 25+1 Lj and color index 1 or 8). (On') are ^ ie l° n g distance nonperturbative matrix 
elements demonstrating the probability of a cc pair evolving into the physical state H. da 
can be obtained from perturbative calculations, while (On') are nonperturbative parameters 
which can not be calculated perturbatively. 

{On / cons i s t °f two kinds of matrix elements, he., the color-singlet and color-octet 
matrix elements. Color-singlet matrix elements may be related to the quarkonium radial 
wave function or its derivatives at the origin, and may be calculated by potential models 
or estimated by leptonic decay widths of quarkonium states. Whereas, the nonperturbative 
color-octet matrix elements can only be determined from lattice QCD calculations or by 
fitting the theoretical prediction of quarkonium production rates to the experimental data. 
The preliminary lattice QCD calculations of quarkonium decay matrix elements have been 
presented in the literature []3J. For the quarkonium production matrix elements, before 
lattice QCD produces its result, the NRQCD velocity scaling rules may be used to give a 
rough (order of magnitude) estimate about size of the color-octet matrix elements. Under 
this velocity scaling rules, the matrix elements (Pn') (including both color-singlet and color- 
octet) are related to each other by orders of m,QV 2 or v 2 , where tuq is the mass of the heavy 
quark and v is the typical relative velocity of the heavy quark in the bound state. For 
example, in J/ip production, the most important NRQCD matrix elements are (of ( 3 S'i)^, 
(p% ( 3 5'i)^, ^(9g’( 1 <So)^ and (p% ( 3 Pj)^), and according to the scaling rules, their relative sizes 
satisfy the following relations 

(pi^Spj ~ m 3 c v 3 , (C$( 3 Si)) ~ m\v\ 

(otCSo)) ~ mlv 7 , (Ot( 3 Pj )} ~ my. (2) 
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Practically, the color-octet matrix elements have been determined by fitting the theo¬ 
retical prediction of quarkonium production to the experimental data at various colliders. 
According to the NRQCD factorization formalism, the matrix elements (On') are universal, 
so that their values measured from different colliders must be the same. This is an impor¬ 
tant test to the color-octet production mechanism. In previous studies |fj] [I|, the matrix 
element is determined by fitting to the high Pt prompt J/ip production at Fer- 

milab Tevatron, and their results are consistent with the velocity scaling rules. However, 
the matrix elements (0^{}Sq)^ and ^(Rg( 3 P 0 )) are not determined individually, but as a 
linear combination ||. Another linear combination of these two elements is also obtained 
in the studies of J/ip photoproduction in e~p collisions [|J and hadroproduction at fix tar¬ 
get experiments [0, and their results are not compatible with that of ||. Recently, some 
further investigations on prompt J/-0 production at the Tevatron are performed || |J. 
Their results show that the parton distributions strongly affect the extraction of the values 
of the color-octet matrix elements, i.e., different sets of parton distributions result in differ¬ 
ent values of the matrix elements. They also show that the initial- and final-state radiation 
are important in prompt J/ij) production at the Tevatron [9;]. All those progresses show 
that the inconsistency between the values measured at CDF and those measured in other 
processes may be mostly due to the large theoretical uncertainties in those calculations, such 
as the parton distributions. 

In contrast, the mechanism of J/-0 production in e + e _ annihilation process is much clearer 
than those hadronic processes discussed above. The parton structure is simpler, and there 
is no higher twist effects to be considered, so the theoretical uncertainty is much smaller. 
Because in e + e _ processes J/^ production has much smaller theoretical uncertainty than in 
hadronic J/ip production processes, it can be used to extract the color-octet matrix elements. 
J/V’ production at e + e _ colliders has been investigated by several authors in literatures [10- 
17]. Braaten and Chen have noted that a clean signature of color-octet mechanism may be 
observed in the angular distribution of J/^j production near the endpoint region at e + e~ 
collider in the low energy region such as at CLEO |16|. In this paper, we make use of the 
previous results of the calculations on the J/^j production in e + e _ process to extract the 
color-octet matrix elements by comparing with the preliminary experimental data at low 
energies [fTH] [|H|. 
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In our previous studies |17|], we have calculated color-singlet and color-octet contributions 
to prompt J/ip production in e + e~ annihilation at different energy scales under the NRQCD 
factorization formalism and we have used the nonperturbative matrix elements which are 
consistent with the NRQCD velocity scaling rules. We find that the color-octet contributions 
dominate over the color-singlet contributions in all energy regions. At low energies, as 
pointed out in |TIJ [T7], the dominant production channels are e + e _ —> J/ijj + g via color- 
octet l S 0 and 3 Pj processes, while at high energies (say, above 20 GeV) the dominant channel 
is e + e~ —> qqg* followed by g* fragmentation to color-octet J/ij) process. 

For the color-octet processes 


e + e -*• 7 * -> g + cc[ 8 , +1 Lj \, 


(3) 


where 2S+l Lj represent the cc pair states 1 Sq and 3 Pj, from Ref. |L 6 |, we readily have 


cr(e + e“ - J/^g) = C s ( OtCS 0 )} + (ot{ 3 P 0 )) , 


(4) 


with 


_ 647r 2 e 2 a! 2 a s 1 — r 

CL = -#- 


s 2 m 


n 2567r 2 e 2 a 2 a s r (l — 3r) 2 | 6(1 + r) | 2(1 + 3r + 6r 2 ). 

Cn — -—o-1—:-1-:-1— 


(5) 

( 6 ) 


9 s 2 m 3 1 — r ' 1 — r ’ 1 — r 

where r = m 2 /s, m is the mass of J/ip, and s is the e + e~ collision c.m. energy squared. 
Here we have used the approximate heavy quark spin symmetry relations 


Oi^Pj)) « (2./ + 1) (ot( 3 P„)) . 


(7) 


To extract the color-octet matrix elements, one needs to subtract the color-singlet con¬ 
tributions from the total cross section. The leading order color-singlet J/-0 production rates 
at low energies (< 25 GeV) comes from the process 


e + e —> 7 * —> J/t/’ + gg- 


( 8 ) 


The cross section of this process is 


da(e + e -> J/'ipgg) 
a^dzdx i 


64e 2 a 2 {pt ( 3 *S'i)) 2 
--- §—~r f(z,xi,r), 


27 


m° 


(9) 
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where 


f{z,xi,r) = 


(2 + x 2 )x 2 
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(2 + xi)x\ 
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( 10 ) 


and 


— <XQED(e + e —> / i + p ). 

The variables z, x, are defined as 


z = 


2 p ■ k 


Xi = 


2 Pi ■ k 


( 11 ) 

s s 

where k, p and pi are the momenta of the virtual photon 7 *, J/^, and the outgoing gluon, 
respectively. 

In the numerical calculations we use the following parameters as input 


m c = 1.5 GeV, a,(2m c ) = 0.26, (of( 3 S,)) = 1.08Gel2 3 [3], 


( 12 ) 


In Fig.l, we show the theoretical prediction of direct J/'ij) production compared with the 
experimental data (note that contribution from ^' feeddown has been subtracted from the 
data sample). The data at yfs = 4 ~ 5GeV are from PLUTO [|H|, and those at yfs = 
10 ,6GeV are from CLEOII [|1^] [|I3| . The two color-octet matrix elements ( 1 S'o)) an d 
( 0 % ( 3 Po)) ar e treated as free parameters to be determined by fitting the data. We use the 
data from two energy points to determine these two matrix elements. One is from PLUTO 
at a/s = 5 GeV and the other from CLEO at yfs = 10.6GeU. At the collision energies 
below 5 GeV there might exist some cc resonances which could contaminate the prompt J/ip 
production, that is why we neglect the data in the energy region below yfs = 5 GeV to carry 
out the extraction. And the extraction result is 


(ptfSf,) = 1.1 x If )~ 2 GeV 3 , 
° t(3P f = 0.74 x 10~ 2 GeV 3 , 


rnt 


(13) 


where the central values of experimantal data (Rf p = 5.0 x 10 3 at fs = 5.0GeU and 
R e f p = 1.0 x 10 ~ 3 at yfs = 10.6GeU) are used and the experimental errors are not included. 
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The dotted line in Fig.l represents the contribution from color-octet 1 S'o subprocess, the 
dashed line is from color-octet 3 Pj process, and their sum is plotted as the dotted-dashed 
line, and the color-singlet contributions is shown as the short dashed line. ^Frorri this figure, 
we can see clearly that the color-singlet contributions alone can not explain the observed cross 
sections in these energy regions. After including the contributions from both color-singlet and 
color-octet production processes (shown as the solid line), a satisfactory agreement between 
theoretical prediction and the experimental data will be achieved. Especially, the cross 
section of direct J/-0 production at y/s = 5 GeV is much larger than that at y/s = 10.6 GeV. 
This behavior can not be explained by the color-singlet model with perturbative QCD. 

Here, we would like to emphasize that the PLUTO data are due to direct J/^j produc¬ 
tion but not the decays of cc resonances Jl8|. Furthermore, we assume that in the energy 
region around y/s = 5 GeV, perturbative QCD is approximately applicable to describe the 
J//> production processes. In this connection, we note that Driesen et al. have applied 
perturbative QCD to calculate the ,J/tjj production in e + e _ process within the color-singlet 
model JT3|. They could explain the experimental data from PLUTO by using a large run¬ 
ning coupling constant, i.e., setting a s (Q 2 ) = 1 if yfQ 1 < 1 GeV. Their results indicate that 
the main contribution to J/rj) production in this process comes from the region where the 
conjunctioned gluons are soft. It is not entirely clear whether this color-singlet cc plus a 
soft gluon can be factorized into a higher Fock state for the physical J/ip, such as color- 
octet 3 Pj states and color-octet 1 S'o state. Nevertheless, if the experimental results of these 
two collaborations (CLEO and PLUTO) are further confirmed, the color-octet production 
mechanism will probably provide a quite unique explanation. 

In the following we estimate the theoretical uncertainties induced by the choice of 
the charm quark mass, the strong coupling constant and the color-singlet matrix element 
Of( 3 Si)\. First, we consider a variation of the charm quark mass m c from lAGeV to 


l.QGeV with other parameters unchanged (as in ([T 2 |) ). This results in a variation of the two 
color-octet matrix elements fitted values, 


(ot^So)) = -0.25 to 2.3 x 10- 2 GeU 3 , 


OtfPo] 


= 0.91 to 0.61 x 10~ 2 GeV 


mi 


( 14 ) 


Also, a variation of the strong coupling constant a s from 0.24 to 0.30 will result in 
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OtPSn)) — 1.8 to 


0 

8 V Po) 


Q^( 3 


0.009 x 10 ~ 2 GeV 3 , 
= 0.76 to 0.71 x 10 ~ 2 GeV 3 . 


m- 


And a variation of the color-singlet matrix element (^Oi( 3 Si] 
will result in 


( 15 ) 

from l.OGeV 3 to 1.2 GeV 3 


(otCSo)) = 1.4 to 0.77 x 10~ 2 GeV 3 , 


Ot( 3 P n 


= 0.72 to 0.76 x 10 ~ 2 GeV 


m- 


( 16 ) 


^Frorn the above results, we can see that the fitted value of the element (0% ( 1 S'o)y is sensitive 
to a variation of the three parameters m c , a s and (O^^Sq)^. This is because, the cross 
section of J/^ production is not sensitive to the value of the element (O^^Sq)^ (he., the 
coefficient of the (O^^Sq)^ term in the expression of cross section is very small, see (f|) 
(H)). In contrast, the fitted value of (of ( 3 Po)) /m 2 remains rather stable when the three 
parameters vary in reasonable ranges. The above rough estimates may show the value of the 
element (0% ( 3 Po)) /K = (0-6 — 0.9) x 10 ~ 2 GeV 3 . 

Finally, we discuss how the present experimental errors can affect the fitted values 
of the two matrix elements. Obviously, the large errors in experimental data will re¬ 
sult in a larger error for the fitted values of the two color-octet matrix elements. To 
see this, in Fig.2 we show four curves for the J/-0 production cross section, where the 
input parameters (a s , m c , (ot{ 3 Si)}) are same as in Fig.l. The two dotted lines cor¬ 
respond to (0% ( 3 Po)) /m 2 c = 7.4 x 10 ~ 3 GeV 3 , (of^So)^ = 5.0 x lCr 2 ChW 3 (up), and 
Og ( 1 A 0 )) = 1.0 x 10 - 4 GeO 3 (down) respectively. The two dashed lines correspond to 
Of^So)) = 1.0 x 10 ~ 2 GeV 3 , (og( 3 P 0 )) /m 2 c = 2.0 x 10- 2 Gef/ 3 (up), and (C^( 3 P 0 )) IK = 
2.0 x 10~ 3 GeI/ 3 (down) respectively. We can see that the total cross section is not sensitive 
to the value of (0% ( 1 5' 0 )^. When it increases from 1.0 x 10 ~ A GeV 3 to 5.0 x 10 - 2 GeId 3 , the 
total cross section only changes a little. But the value of (0% ( 3 Po)) varies rapidly as the 
data value changes. Therefore, in order to determine the color-octet matrix elements more 
rigorously, more precise measurement for the J /-0 cross section at low energies is apparently 
needed. 

Meanwhile, we note that the color-octet matrix elements discussed in this paper have 
also been determined from other experiments. At the Tevatron, prompt J/-0 production 
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at high Pt has been compared with theoretical predictions, and a global fit to all Pt region 
shows that at low Pt the theoretical prediction is dominated by the contributions from 
< Of ( 1 *S'o) > and < Of( 3 Po) >, and the fitted result is J5J 


o* (‘So)) + P (O 8 *( 3 fo)) = 6.6 x 10~*GeV 3 . 


( 17 ) 


Recently, Beneke et al. || updated the extraction of color-octet matrix elements by com¬ 
paring the Pt distribution for unpolarized direct J/'ijj production with the most recent CDF 
data 0- By using the CTEQ4L PDF (parton distribution function) set (in their paper, 
they have used three PDF sets, which resulted in different extractions ||), they obtained a 
different linear combination of the two color-octet matrix elements, 


(ofCSo)) + ^ (ot( 3 Po)) = 4.38 x 10 - 2 GeR 3 . ( 18 ) 

m c 

In Ref. 0, Sanchis-Lozano et al. used a Monte Carlo event generator to treat the high 
order initial- and final-states radiation. By fitting to the experimental data of prompt J/-0 
production at the Tevatron, the authors obtained another extraction of these two color- 
octet matrix elements, 


OtCSo)) + -4 M( 3 P 0 )) = 1.44 x 10- 2 GeR 3 , 

/ 777 “ \ / 


(19) 


where the CTEQ2L PDF set is used. The value in Eq.(19) is much smaller than those in 
Eq.(pT|) and (|18D. This shows that the inital-states radiation effects are crucialy important 
to the extraction of these two color-octet matrix elements at the Tevatron. 

The studies of photoproduction at e~p collisions show that the matrix element values in 
Eq.(|T7|) may be overestimated [|(|, and the authors obtained another linear combination of 
these two elements from the forward J/ijj photoproduction cross-section measurements, 


(ofCSo)) + ^ (ot( 3 P 0 )) = 2.0 x 10- 2 GeE 3 . (20) 

The fixed-target hadroproduction result [P] gives the same argument against the matrix 
elements values in Eq. ([17]) , and gives 


OtCS 0 )) + = 3.0 x 10- 2 GW 3 . 




( 21 ) 


Our results (|T^) are not compatible with (^0|) and (^T[). The inconsistency may be mostly 
due to theoretical uncertainties and large experimental errors. 
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However, it should be noted that there are large theoretical uncertainties in our extraction 
of the color-octet matrix element ( 0 % ( 1 S'o)^), and also the large experimental errors may 
further affect the extraction of the two matrix elements ^(9g( 1 S'o)^ and <^Og( 3 Po))- So, 
the compatibility of our results with those from Tevatron and the incompatibility of our 
results with those from fixed-target hadroproduction and photoproduction should not be 
taken seriously. 

In Ref. B Fleming et al. make an analysis of the present extraction values of the 
color-octet matrix elements. After considering the constraints by the J/ip production rate 
in 6-decays, they obtained a rough range of the linear combination of the two color-octet 
matrix elements, 


[ofCSo)) + — (ofCPo)) = 1.0 ~ 6.0 x 10~ 2 GeV 3 . 

m c x 

This is compatible with all the above results (Eqs.([L3]) and ([171) through (|2l|) ). 


( 22 ) 


When this work is in progress, we receive a preprint of Fleming et al. |22|], in which 
they calculate the J/if) leptoproduction. They find that a negative value for the color-octet 
matrix element (0% ( 3 Po)) is still possible, which is inconsistent with our results. We would 
like to point out that if this matrix element is negative, the color-octet 3 Pj processes would 
provide a large negative contribution to J/'ijj production in e + e _ process at low energies due 
to a large positive coefficient for the matrix element (0% ( 3 Po)) a t l° w energies (see Eqs.(4)- 
(6)). If so, the experimental data in Fig.l can not be explained. The conflict between our 
results with those of |Z| about the color-octet matrix elements puts forward a question to 
the universality of the color-octet matrix elements in these two situations. 

In conclusion, in this paper we have calculated the direct J/^j production rates in e + e~ 
process at low energies. Our results show that the color-octet production mechanism is 
significant to the production cross section in this energy region, and this may provide another 
positive test for the validity of the color-octet production mechanism. The extraction of 
elements (0% ( 1 *S'o)^ and ( 3 Pj)^ are performed by fitting to the experimental data. But 
the large errors of the preliminary experimental data can not give a precise estimate of the 
color-octet matrix elements. We hope that more precise experimental data of direct J/ip 
production in e + e _ annihilation at CLEO II and BEPC will soon be obtained. 
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Figure Captions 


Fig.l. The total cross section of J/-0 in e + e^ process. Color-octet 1 S , 0 process contribution 
is represented by the dotted line, color-octet 3 Pj process by the dashed line, the sum of 
these two color-octet processes by the dotted-dashed line, and color-singlet contribution by 
the short dashed line. The experimental data are taken from Ref. |1{J Jl8| . 


Fig.2. Description of the allowed values of the color-octet matrix elements constrained by 
the experimental data. The dotted lines correspond to (0% ( 3 Po)) / m c = 7A x 10 -3 GeF 3 , 
(0%(}Sq)^ = 5.0 x 10~ 2 GeF 3 (up), and (0%(}Sq)^ = 10 -4 GeF 3 (down) respectively. 
The dash lines correspond to (0%( 1 Sq)^ = 1.0 x 10~ 2 GeF 3 , (0% ( 3 Po)) /rn 2 c = 2.0 x 
10” 2 GeF 3 (up), and (Cg ( 3 P 0 )) K = 2.0 x 10 -3 GeF 3 (down) respectively. 


12 





^(XIO- 4 ) 



Fig-1 





^(XIO- 4 ) 




